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ABSTRACT 

The magnetic field of the solar corona has a large-scale dipole character, which maps into 
the bipolar field in the solar wind. Using standard representations of the coronal field, we show 
that high-energy ions can be trapped stably in these large-scale closed fields. The drift shells 
that describe the conservation of the third adiabatic invariant may have complicated geometries. 
Particles trapped in these zones would resemble the Van Allen belts and could have detectable 
consequences. We discuss potential sources of trapped particles. 

Subject headings: acceleration of particles — Sun: corona — Sun: magnetic fields 



1. Introduction 

The large-scale magnetic field of the solar 
corona may trap high-energy particles for long 
perio ds of time, according to early suggestions 
(e.g.. lElliotl [l969h . The original impetus for con- 
sidering coronal trapped particles came from the 
presence of long timescales (hours to days) and 
the appearance of diffusive transport of solar en- 
ergetic particles (SEPs, then often called "solar 
cosmic rays"). Diffusive transport suggested coro- 
nal trapping and "leakage' into the heliosphere. 
More recently it has become clear that SEPs 
can be accelerated in heliospher ic shock fronts 
(|Cane et al.l Il988t iReamed [l999h . and that the 
long timescale for injection into the heliosphere 
could be associated with the propagation time of 
the shock. Because the shock disturbance takes 
hours to days en route to 1 A.U., a leaky coronal 
trap now seems unnecessary to explain SEPs. It 



should be noted, however, that observations of the 
2.223 MeV 7-ray line due to neutron capture have 
already suggested trappin g for many hours even 
in active-region fields (e.g. Ryan 200(1) . 

There are several other implications of coronal 
particle storage, to the extent that p lasma insta- 
bilities permit it (e.g.. IWentzel 1976). In partic- 
ular, as noted by Elliot and others, the particle 
storage may involve a substantial energy density 
and perpendicular pressure. The particles could 
therefore have effects on the structure of the coro- 
nal plasma. However the ion component, not cou- 
pling well to detectable electromagnetic radiation, 
would mainly be identifiable via its secondary ef- 
fects. 

The geometry of the large-scale coronal mag- 
netic field, discussed below in Section [2 strongly 
suggests a basically dipole character at solar mini- 
mum and in the high corona. If this were true, 



then we would expect the existence of "forbid- 
den zones" in St0rmer's sense, with the possi- 
bility of long-term trapping of particles result- 
ing from cosmic or solar sources. In this Letter 
we revisit the geometry of coronal particle trap- 
ping and demonstrate that it can be complete, 
in the sense that all three of the adi abatic in- 
varia nts of charged-particle motion (e.g. Northropl 
19631 ) may be conserved for times long enough 
for particles to circulate completely around the 
Sun and thus potentially complete closed trap- 
ping shells. Section [2] describes the methods used 
to establish this result, which is based on anal- 
ysis and modeling. The deta iled modeling uses 
the lSchriiver fe DeRosa (2003) version of the "po- 
tential field source surface " (PFSS) approach to 
the solar large-scale fields ([Altschuler fc Newkirk 
19691 : ISchatten et"aD Il969h . This type of model 
uses the routine Zeeman-splitting observations of 
the line-of-sight photospheric magnetic field as an 
inner boundary condition, and assumes a radial 
field outside a fixed "source surface" at (in this 
case) 2.5 R©. In between the photosphere and the 
source surface, the model simply uses a potential 
expansion of the field; the body currents known 
to flow in the corona are replaced by currents out- 
side the source surface. This approach appears to 
work reasonably well for a large-scale structure. 
In Section [3) we summarize the experiments we 
have made using PFSS fields. Even in these rel- 
atively complicated fields we find trapping orbits 
that allow particles to circulate all the way round 
the Sun, pote ntially conservin g the third adiabatic 
invariant (e.g. Northropl 19631 ). 



2. The coronal magnetic field 

Eclipse pictures show the solar-minimum corona 
to have a strongly dipolar configuration, especially 
at solar minimum. Hale's discovery of solar mag- 
netism allowed an explanation of this appearance 
by analogy with the Earth's field, which comes 
from sources deep in the interior and has a true 
dipole nature. The coronal field has a fundamen- 
tal difference in that at low altitudes small-scale 
fields dominate; these fields originate in the mag- 
netic active regions or in the "net work" structure 
(e.g., IZwaanl [l987t lHarvevI fl993h . The dipolar 
character only appears at higher altitudes in the 
corona, where the intense small-scale fields of solar 
active regions and the network do not extend. The 



dipole structure thus appears only on the largest 
scales, and is mainly visible in eclipse images at so- 
lar minimum. The observed polarity structure of 
the magnetic field in the so lar wind normally ha s 
a simple "sector" structure ( Wilcox fc Ness 19651) . 
It is thus bipolar although highly non-potential in 
character since the field is largely radial. 

Our detailed knowledge of the coronal field re- 
mains meager, with almost negligible direct mea- 
surement because of the difficulty of making pre- 
cise polar i zation measurements in the corona (e.g. 



Lin et al. 2000| ). There are several approaches 



to extrapolating the photospheric field into the 
corona instead. These can be compared with in 
situ measurements of the heliospheric field, but 
have the drawback that the treatment of cur- 
rents penetrating the photosphere and flowing 
in the corona is ambiguous even in the "non- 
linear fojxe^^e_^gJd"_(NLFFF) extrapolations 



(e.g., ISchriiver et alj [2006h . Full MHD simula 



tions of coronal structure (e.g.. iMikic et al. 2007 ) 
may improve our detailed understanding. In 
the meanwhile, the Schrijver & DeRosa PFSS 
computational structure is very convenient, be- 
cause it runs as a softwa re package in SolarSoft 
(jFreeland fc Handvlll9981 ) and the database gives 
6-hour time resolution since 1996, based on Solar 
and Heli ospheric Observatory /Michelson Doppler 



ry/. 

Imager ([Scherrer et al.l 119951 ) . We can use this 



model to explore particle trapping in the coronal 
field, but note that it has physical limitations. 
Indeed, neither PFSS nor more any of the more 
sophisticated coronal field models have succeeded 
in reliably detecting the obvious restructurings 
due to flares and coronal mass ejections (CMEs), 
which have clear pho tospheric counterparts (cf. 
Sudol fc Harvevll2005f) . 



3. Particle propagation 

In a low-order analytic model of the field such 
as a dipole or a Mead fielcfl, the length scales of the 
field greatly exceed the particle gyroradii, even for 
the heaviest ions we consider. The same should be 
true of the PFSS representation or other extrapo- 
lations of the photospheric field with limited spa- 



1 The Mead fields llMeadlll964l) are distorted dipoles with a 
full analytical description, originally intended as a repre- 
sentation of geomagnetic field. We use these to test our 
codes. 
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tial sampling. We therefore solve the relativistic 
gu iding-center drift eq uation s of m otion as given 
bv lNorthropI (|l963l ) or lParkd (|2004h . using spher- 
ical polar coordinates. From the point of view 
of particle trapping in the corona, the question 
is whether or not particles can readily move via 
gradient or curvature drifts from closed to open 
field lines, and hence "escape," and whether their 
mirror points move low enough for precipitation. 
Particles that can survive either of these threats or 
collisional energy decay then can drift completely 
around the Sun, potentially conserving the third 
adiabatic invariant and defining a drift shell. 

Our modeling here aims at defining the mor- 
phology of the particle trapping, rather than pro- 
viding quantitative results. We follow particle 
motions with an RK4 (Runge-Kutta) integrator 
(jPress et al.lll986j ) as implemented in IDL, for con- 
venience, and typically use a heavy test parti- 
cle (mass number 207, for example) to exagger- 
ate the cross-field drifts. The collisional energy 
losses in the corona l plasm a are estimated from 
Weaver fc Westprial (2002), who provide correc- 



tions to the basic formula —dE/dx ex Z 2 n e L/v 2 , 
where Z is the projectile charge, v is its velocity, 
n e is the electron density, and L is the "stopping 
logarithm." For an ionized medium we take the 
ionization potential to be hujp, where h is Planck's 
constant, and uj p is the plasma frequency. We 
adopt the Van de Hulst-Newkirk density model 
(jNewkirkl 119611 ). and Figure Q] shows representa- 
tive lifetimes for 1-100 MeV protons. The coronal 
magnetic field is not really understood well enough 
to justify detailed calculations; the main concern 
here is that the PFSS approach, though conve- 
nient, systematically misrepresents it both in the 
low c orona where field-alig ned currents are strong 
(e.g., ISchriiver et al. 2008) and also in the high 
corona near the ad hoc source surface at 2.5 R . 
Thus we would trust our numerical experiments on 
trapping mostly in the middle corona, at heights 
of 1.5-2.0 R©. It would clearly be possible to do 
similar experiments with more complete models, 
such as those bas ed on the MHD equations (e.g. 
Mikic et al.1l2007t ). 



4. Existence and stability of trapping 
zones 

After testing our code on analytical field mod- 
els (dipole and Mead) , we have explored parameter 
space by injecting particles with specified energies 
and pitch angles into PFSS coronal models. The 
initial experiments were with the PFSS field real- 
ization of 2006 September 30 12:04 UT, selected 
just for its simple appearance (see Figure [5]). The 
easily seen non-dipolar structure shows that this 
epoch is not the very simplest that could have been 
chosen, and indeed four GOES X-class flares oc- 
curred subsequently in 2006 December. 

Figure [3] shows the trajectory of a single test 
particle injected at a chosen point (near the mid- 
dle of the estimated trapping zone in this case). 
We used an unusual ion (Bi +83 ) to exaggerate the 
losses. The bouncing motion results from the fact 
that the injection point (height 2.0 Rq) was at 
the apex of the field line intersecting the injec- 
tion point, rather than its minimum value of |B|. 
This particle completed its full shell motion in 
about 6 hr and could be followed for more than 
one circulation. It could be argued that the PFSS 
model is not realistic enough, simply because the 
field reconstruction is limited to a small number of 
spherical harmonics, and hence cannot show steep 
gradients and large non-adiabatic effects. On the 
other hand we are using a test ion with a Lar- 
mor radius orders of magnitude too large, and we 
are experimenting around 2 Rq where the field 
presumably cannot be stressed strongly on small 
scales. We thus feel that this approximation is re- 
alistic enough to establish the principle that trap- 
ping zones exist in the corona. 

How large is the volume of phase space show- 
ing this property? Figure |4] shows the results for 
test particles in a half of the meridional plane (9, 
R-sp&ce) at the injection longitude (6.02 rad) for 
the particle of Figure [31 Particles injected at zero 
pitch angle (but not necessarily zero parallel veloc- 
ity) have trapped orbits, and the remainder either 
start on open field lines or drift onto them, with 
a few particles at the lower boundary precipitat- 
ing. We conclude from this and other tests that 
drift shells (stable trapping) occupy a substantial 
fraction of the coronal phase space for field lines 
closing above 1.4 R Q . The drift-shell regions may 
not be contiguous. 
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Fig. 1. — Lifetimes (E/E) of energetic protons in 
a representative PFSS field mode( 2006 Sep. 30)1, 
for different injection altitudes at (9, <fr) = (6.02, 
1.04). The sets of points correspond to proton 
energies of 4, 8, 16... 256 MeV. 




Fig. 3. — Trajectory of an individual test particle 
(Bi +83 to exaggerate the non-adiabatic effects) in 
the PFSS field of Figure H The initial energy is 
150 MeV, pitch angle 90°, and coordinates 0, 
R) = (6.02, 1.04, and 2.0 R ). 




Fig. 2.— PFSS field for 2006 September 30, 12:04 
UT. Open field lines are in color, closed in white. 




Fig. 4. — Sector of the meridional plane 
(60° x 1.5 R©) at heliolongitude 6.02 rad for the 
same date as the simulations shown above. The 
plus signs mark the long-lived orbits for particles 
injected near the tops of the field lines intersecting 
the grid of points, with zero initial parallel veloc- 
ity. 
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5. Sources of particles 

In the Earth's radiation belts, the "Cosmic 
Ray Albedo Neutron Decay" (CRAND) mecha- 
nism provides a basic guarantee that some par- 
ticles can materialize in the Van Allen trapping 
zones. The same is true for the Sun, except that 
the cosmic-ray fl ux is diminished b y heliosphcric 
modulation (e.g. ISeckel et al.l Il99lh - In addition 



the mean atomic number of the solar atmosphere 
is lower, around 1.3 versus 7.2 for the Earth's at- 
mosphere. This means more (p,p) primary inter- 
actions for the primary cosmic rays, which pro- 
duce fewer neutrons than do spallation reactions 
on higher- Z nuclei. 

There are good reasons to think that trapping 
zones in the coronal magnetic field will actually 
contain particles. Large-scale shock waves propa- 
gate through the region as a resu lt of flares and /or 
CMEs. Shock acceleration (e.g. Reames 1999f) is 
known to be an important factor in the produc- 
tion of SEPs. We cannot be quantitative now 
about the possible contributions of this mecha- 
nism to the trapping zones, especially since the 
CME closely involved with causing the shock, 
and the particle acceleration, itself may disrupt 
the geometry of the corona. Charge exchange 
at MeV energies recently has been observed via 
flare-associ ate energetic neutral a toms observed by 
STEREO (jMewaldt et al.l l2009h . Although only 
a single event (2006 December 5) has thus far 
been reported, this clear observation strongly sug- 
gests that SEP-associated flares will populate dis- 
tant regions of the corona via neutralization and 
reionization interactions. Similarly, the copious 
production of neutrons by flares, as evidenced by 
the 2.223 MeV 7-ray emission and by their di- 
rect detection at Earth, guarantees a source of 
particles at some level. Neutron-decay particles 
have also been detect ed directly in the heliosphere 
( Evenson et al. I ll983h . Other mechanisms operate 



in the Earth's radiation belts, and via pitch-angle 
diffusion the trapped particles can actually gain 
energy as they migrate to regions of stronger field. 

6. Conclusions 

We have made basic models of high-energy par- 
ticle propagation in the solar corona, using PFSS 
extrapolations of the magnetic features observed 
at the photospheric level. The PFSS models, while 



not perfect, have a reasonable track record for es- 
tablishing connectivity on large scales. Our cal- 
culations show that trapping zones exist in the 
corona, just as they do in the Earth's field. In 
these zones the main "escape route" for a trapped 
particle will be its collisional losses at lower ener- 
gies. 

While it is not possible now to study these 
particles in situ, as one can do near the Earth, 
their presence may have indirect signatures. 
These could include 7-ray, X-ray, or radio emis- 
sions, and the particles could also have dy- 
namical consequences for the structure of the 
plasma if numerous enough. Direct remote sens- 
ing of the ion component is almost nonexis- 
tent, but high-energy electrons could possibly 
be detectable. Variou s mechanisms exi s t (e.g . 
Moskalenko et all l2006t iKahler fc Ragotl l2008h . 
and both EGRET (jOrlando k Strong) l2008h as- 
yet-unpublished Fermi results show a distinct 
source of high-energy 7 rays from the solar di- 
rection. The particle environment of the corona 
itself may play a role in understanding this source 
and its relationship to solar-system background 
sources for Fermi observations. 
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